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A STUDY OF THE MOTION AND AERODYNAMIC HEATING OF BALLISTIC MISSILES ENTERING 
THE EARTH’S ATMOSPHERE AT HIGH SUPERSONIC SPEEDS 1 

By H, Julian Allen and A. J. Eggers, Jr. 


SUMMARY 

A simplified analysis is made of the velocity and deceleration 
history of ballistic missiles entering the earth's atmosphere at 
high supersonic speeds. It is found that , in general } the gravity 
force is negligible compared to the aerodynamic drag force and, 
hence , that the trajectory is essentially a straight line. A con- 
stant drag coefficient and an exponential variation of density with 
altitude are assumed and generalized curves for the variation of 
missile speed and deceleration with altitude are obtained. A 
curious finding is that the maximum deceleration is independent 
of physical characteristics of a missile {e. g ., mass, size, and 
drag coefficient) and is determined only by entry speed and 
flight-path angle, provided this deceleration occurs before impact. 

The results of the motion analysis are employed to determine 
means available to the designer for minimizing aerodynamic 
heating. Emphasis is placed upon the convective-heating 
problem including not only the total heat transfer but also the 
maximum average and local rates of heat transfer per unit area. 
It is found that if a missile is so heavy as to be retarded only 
slightly by aerodynamic drag , irrespective of the magnitude of the 
drag force, then convective heating is minimized by minimizing 
the total shear force acting on the body. This condition is 
achieved by employing shapes with a low pressure drag. On the 
other hand, if a missile is so light as to be decelerated to rela- 
tively low speeds , even if acted upon by low drag forces, then 
convective heating is minimized by employing shapes with a 
high pressure drag, thereby maximizing the amount of heat 
delivered to the atmosphere and minimizing the amount delivered 
to the body in the deceleration process. Blunt shapes appear 
superior to slender shapes from the standpoint of having lower 
maximum convective heat-transfer rates in the region of the nose. 
The maximum average heat-transfer rate per unit area can be 
reduced by employing either slender or blunt shapes rather than 
shapes of intermediate slenderness. Generally, the blunt shape 
with high pressure drag would appear to offier considerable 
promise of minimizing the heat transfer to missiles of the sizes, 
weights, and speeds of usual interest. 

INTRODUCTION 

For long-range ballistic trajectories one of the most diffi- 
cult phases of flight the designer must cope with is the re- 
entry into the earth's atmosphere, wherein the aerodynamic 
heating associated with the high flight speeds is intense. 
The air temperature in the boundary layer may reach values 
in the tens of thousands of degrees Fahrenheit which, com- 

• Supersedes NACA Technical Note 4047 by H. Julian Allen and A. J. Eggers, Jr. ( 1957. 


bined with the high surface shear, promotes very great 
convective heat transfer to the surface. Heat-absorbent 
material must therefore be provided to prevent destruction 
of the essential elements of the missile. It is characteristic 
of long-range rockets that for every pound of material which 
is carried to “burn-out," many pounds of fuel are required 
in the booster to obtain the flight range. It is clear, there- 
fore, that the amount of material added for protection from 
excessive aerodynamic heating must be minimized in order 
to keep the take-off weight to a practicable value. The 
importance of reducing the heat transferred to the missile 
to the least amount is thus evident. 

For missiles designed to absorb the heat within the solid 
surface of the missile shell, a factor which may be important, 
in addition to the total amount of heat transferred, is the 
rate at which it is transferred since there is a maximum 
rate at which the surface material can safely conduct the 
heat within itself. An excessively high time rate of heat 
input may promote such large temperature differences as 
to cause spalling of the surface, and thus result in loss of 
valuable heat-absorbent material, or even structural failure 
as a result of stresses induced by the temperature gradients. 

For missiles designed to absorb the heat with liquid 
coolants (e. g., by transpiration cooling where the surface 
heat-transfer rate is high, or by circulating liquid coolants 
within the shell where the surface heat-transfer rate is 
lower), the time rate of heat transfer is similarly of interest 
since it determines the required liquid pumping rate. 

These heating problems, of course, have been given con- 
siderable study in connection with the design of particular 
missiles, but these studies are very detailed in scope. There 
has been need for a generalized heating analysis intended to 
show in the broad sense the means available for minimizing 
the heating problems. Wagner, reference 1, made a step 
toward satisfying this need by developing a laudably simple 
motion analysis. This analysis was not generalized, how- 
ever, since it was his purpose to study the motion and heat- 
ing of a particular missile. 

It is the purpose of this report to simplify and generalize 
the analysis of the heating problem in order that the salient 
features of this problem will be made clear so that successful 
solutions of the problem will suggest themselves. 

A motion analysis, having the basic character of Wagner's 
approach, precedes the heating analysis. The generalized 
results of this analysis are of considerable interest in them- 
selves and, accordingly, are treated in detail. 
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Equations (1) then reduce to the single equation 

d'y dV , CvPoA -py* 

dT> = ~dt' 9+ 2m 


analysis 

motion of the body 

Consider a body of mass m entering the atmosphere from 
file ***>’ the 



Noting that 


we let 


2m 
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( 3 ) 
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„d equation (3) becomes the linear differential equation 
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As an example, consider the vertical descent o*^ ^ 
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Figure 1. — Variation of density with altitude. 
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Figure 2. — Variations of velocity with altitude for a 1-foot diameter, 
solid iron sphere entering the earth’s atmosphere vertically at 
velocities of 10,000, 20,000, and 30,000 ft/sec. 


vertical entrance velocities of 10,000, 20,000, and 30,000 feet 
per second at 40 miles altitude which, for these cases, may 
be considered the “outer reach” of the atmosphere. The 
results of these calculations are presented as the solid curves 
in figures 2 and 3. 

It is seen in figure 3 that for the high entrance speeds con- 
sidered, the decelerations reach large values compared to 
the acceleration of gravity. This suggests that the gravity 
term in equation (3) may be neglected without seriously 
affecting the results. 5 When this term is neglected the 
equation of motion becomes 
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Figure 3. — Variations of deceleration with altitude for a 1-foot 
diameter, solid iron sphere entering the earth’s atmosphere vertically 
at velocities of 10,000, 20,000, and 30,000 ft/sec. 

Integration gives 

in e ~ fivJ f~ const - 

or 

- c-Sv 

V=const.Y.e Vm 

At the altitude of 40 miles it can readily be shown that the 
term 

_£o^4 

e 

is very nearly unity so that the velocity may be written 

_£dp^A t _ gv 

V=V B e (9) 


1 It Is usual to neglect the gravity acceleration a priori (see e. g., refa. I and 4.) 
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and the value of the maximum deceleration is 

/dV\ /dV 


dV 


W T Cpf> 0 A -fly 

dt CdPoAVe 2 r -&vp 


0m 


( 10 ) 


where V* is the entrance speed. 

Bv use of equations (9) and (10) the vertical-descent speeds 
and^decelerations for the 1 -foot-diametqr sphere previously 
considered have been calculated for tile same entrance > speeds. 
The results are shown as the dashed curves in figures 2 and 3. 
I Ts set that these approximate calculations agree very well 
lilh those based on the more complete equation of motion 

(C The above finding is important, for it indicates that in ^the 
general case, wherein the body enters the atmosphere^ 
high speed at angle B s to the horizontal, the giavitv ter , 
provided 6 E is not too small, may be neglected in equation 
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If equations (13) and (14) are rewritten to make the 
altitude reference point y, rather than zero, then 
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roQTFPctivelv where Ay is the change in altitude from y x . 
sXSltu 5 05) into these expression, can 

readily be shown to give 


' Ay 


and 


so that the flight path is essentially a straight line (i. e., 
e=e E ), and the resultant deceleration equation becomes 

_dy_ C D pAV* (12) 

dt 2m 

Now, again, if the density relation given by equation (2) is 
used and it is noted that 
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Eauations (18) and (19) are generalized expressions foi vdoc- 
itv and deceleration for bodies of constant drag coefficient 
and together with equations (15) and (17), can be used to 
determine the variation of these T'antities wit ^ ™ 

specific cases. The dependence of F (fi&v) and F (fiAy) 
P&y is shown in figure 4. 
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The altitude y, at which the maximum deceleration occurs 
is found from this relation to be 


Figure ^.-Variations of F'tfi Ay), F" 0 Ay), and F'" (fiA V) with 
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If y, is positive the velocity V x (from cqs (13) and (15)) at 
which the maximum deceleration occurs becomes 


The maximum deceleration and the velocity for maximum 
, I“ U s given by equations (17) and (16) apply only 

ff' the altitude y„ given by equation 

wise the maximum deceleration in flight occurs at 

with the velocity (see eq. (13)) 
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HEATING OP THE BODY 


it IK iZ Tfr V1 °r ly IT f ° r practbabla r ocket missiles, 
"rr at e we ? ht 0 the missile be kept to a minimum. 
a V le beat transferred to a missile from the air must be 

■ ; absorbed by s °me "coolant” material. Since this material 

has a maximum allowable temperature, it follows that it 
can accept only a given amount of heat per unit weight 
; Hence, the total heat input to the missile must be kept at a 
minimum for minimum missile weight. 

Often the coolant material is simply the shell of the missile 

■ and as such must provide the structural strength and rigidity 

. for the missile as well. The strength of the structure is dic- 

? !" part > by tbe stresscs induced by temperature gradi- 

■ Z * Wl t hm , . e sh , e ' ■ Slnce thesc romperature gradients are 
proportional to the time rate of heat input, the maximum 
time rate of heat input is important in missile design. The 
hea mg, of course, varies along the surface but, since the 
shell transmits heat along as well as through itself the 
strength of the structure as a whole may be determined by 
the maximum value of the average heat-transfer rate over 
the surface. This is simply the maximum value of the time 
rate of heat input per unit area. On the other hand the 
structural strength at local points on the surface may be de- 
termined primarily by the local rate of heat input. Hence 

e maximum time rate of heat input per unit area at the sur- 
face element where the heat transfer is greatest may also be 
of importance in design. 

If liquid cooling is employed, the maximum surface heat- 
transfer rates retain their significance but, now, in the sense 
ey dictate such requirements as maximum coolant 
pumping rate, or perhaps shell porosity as well in the case of 
t^pintion COOl T g ' Whichever the case > in the analysis 
treated ’ elementS of the heatin g Problem will be 

1. The total heat input. 

2. The maximum time rate of average heat input per unit 
area. 

3. The maximum time rate of local heat input per unit area 
bince it is the primary function of this report to study 

means available to the imssile designer to minimize the heat- 
mg problem, the analysis is simplified to facilitate compari- 
son of the relative heating of one missile with respect to 
■ , be , r accurate determination of the absolute heating of 

mfnd th Ml 69 59 n0t attempted - With this poinf in 
are made followin S assumptions, discussed in A D pendix B, 

J ' Convecti ^ e beat transfer predominates (i. e., radiation 
effects are negligible). 

2. Effects of gaseous imperfections may be neglected. 
glected 10Ck ~ WaV<? boundarydayer interaction may be ne- 

4. Reynolds’ analogy is applicable. 

5. The Prandtl number is unity. 


fer T fr taI “ put -~ Tbe tlmc ra te of convective heat trans- 
fer from the air to any element of surface of the body may be 
expressed by the well-known relation 
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It is convenient in part of this analysis to determine the 
heating as a function of altitude. To this end, noting that 

dt=^L- 

V sin d E 

we see that equation (22) may be written 
dH h t {T r -Tj t 
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V sin 0, 
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With the assumption that the Prandtl number is unity the 
recovery temperature is 

T r=T, (l+l=iM,^=r(l+l=ll il/ J ) 

where 

M Mach number at the altitude y, dimensionless 
y the ratio of specific heat at constant pressure to that 
at constant volume, C,/C„ dimensionless 
1 static temperature at the altitude y, °R 

so that 
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It is seen that for large values of the Mach number, which is 
the case of principal interest, the third term is large com- 
pared to reasonably allowable values of T—T It will 
therefore be assumed that T— T w is negligible • so”that 


Moreover, since 
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equation (24) may be written 
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high speed during descent, the assumption Is obviously remissible EveiTforhlT ' ’ ntal h 
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The total heat input to the body at impact follows from equa- 
tion (29) (integrating over the limits 0<y< 00 ) and is 


6 

Now the local heat-transfer coefficient h t is, by Reynolds’ 
analogy, for the assumed Prandtl number of unity 


hi — 2 


(26) 


where C fl is the local skin-friction coefficient based on con- 
ditions pi, V t , etc., just outside the boundary layer. Thus, 
since (T r —T„) is essentially constant over the entire surface 
S, the rate of total heat transfer with altitude becomes from 
equations (23) through (26) 
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where Q is the heat transferred to the whole surface S. This 
equation may be written 
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(33) 


dQ_ C/pVVS 


dy' 
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wherein C, Vl is set equal to C„ and 

v-aiMM)* (28) 

The parameter C/ is termed "the equivalent friction coeffi- 
cient/ 7 and will be assumed constant, 7 independent of alti- 
tude, again on the premise that relative rather than absolute 
heating is of interest. With equations (2) and (13), then, 
equation (27) is written 


C [)P 0 A . 

<IQ C/ S p 0 V e 1 

dy~ 4sin0 B 


(29) 


Comparison of equation (29) with equation (14) shows that 
the altitude rate of heat transfer is directly proportional to 
the deceleration, so that 


dQ/dy mg / C/S 
( d Vjdt ^ 2 sin \CV4 


) 


(30) 


and therefore the maximum altitude rate of heat transfer 
occurs at the altitude y x (see eq. (15)) and is given by 

/dQ\ $™vj (c;s\ ( l) 

4 —KcoA) } 

It follows, of course, that the altitude rate of heat transfer 
varies with incremental change in altitude from y\ in the same 
manner as deceleration, and thus (see eq. (19)) 


(dQ/dy ) , 


(32) 


i This assumption would appear poor at first glance since the Mach number and Reynolds 
number variations arc so large. Analysis has indicated, however, that the effects of Mach 
number and Reynolds number variation are nearly compensating. The variation in <7/ for 
typical conical missiles was found to be, at most, about 50 percent from the maximum <7/ in 
the altitude range in which 80 percent of the heat Is transferred. 


The impact velocity, V 0 (the velocity of body at y— 0), is 

F„=V*e 2ffmsln, « 

so that equation (33) may be written in the alternative form 


(34) 


Maximum time rate of average heat input per unit area.— 

To determine the time rate of average heat transfer per unit 
area, equations (25), (26), and (28) with equation (22) may 
be shown to give , 

C/pV 3 (35) 

at 

which, together with equations (2) and (13), becomes at 
altitude y 

3 C DP n A —fly 

dH at C/ PqVe* c -BV c 20m sin 9 (36) 

dt 4 

The maximum time rate of average heat transfer per unit area 
is found from this expression to be 


/ dH a A i 
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and it occurs at the altitude 


2 / 2 = In 
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where the velocity is 


v 2 =v E e-*^nv E 


(37) 


(38) 


(39) 


As with altitude rate of heat transfer, it can be shown that 

<4o> 

Equations (37), (38), and (39) apply if the altitude for maxi- 
mum time rate of average beat transfer per unit area occurs 
above sea level. If y 2 , by equation (38), is negative, then 
this rate occurs at sea level and is, from equation (36), 


/ dHg A _/ dHg A _ C/p 0 Ve 3 g - ^msin 
\ dt \ dt J 0 4 


(41) 


Maximum time rate of local heat input per unit area. — 

The elemental surface which is subject to the greatest heat 
transfer per unit area is, except in unusual cases, the tip of 
the missile nose which first meets the air. It seems unlikely 
that a pointed nose will be of practical interest for high-speed 
missiles since not only is the local heat-transfer rate ex- 
ceedingly large in this case, but the capacity for heat reten- 
tion is small. Thus a truly pointed nose would burn away. 
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Body shapes of interest for high-speed missiles would more 
probably, then, be those with nose shapes having nearly 
hemispherical tips. The following analysis applies at such 
tips. 

It is well known that for any truly blunt body, the bow 
shock wave is detached and there exists a stagnation point 
at the nose. Consider conditions at this point and assume 
that the local radius of curvature of the body is a (see sketch). 



The bow shock wave is normal to the stagnation streamline 
and converts the supersonic flow ahead of the shock to a low 
subsonic speed flow at high static temperature downstream 
of the shock. Thus, it is suggested that conditions near the 
stagnation point may be investigated by treating the nose 
section as if it were a segment of a sphere in a subsonic flow 
field. 

The heat-transfer rate per unit area at the stagnation 
point is given by the relation 

NUrkrjTy—Tr) 

(it <T 

where k r is the thermal conductivity of the gas at the re- 
covery temperature (i. e., total temperature) T r , and Nu r is 
the Nusselt number of the flow. If the flow is assumed to 
be laminar and incompressible, 8 Nu r is given, according to 
reference 5, by the relationship 

Nu t = 0. 934 Re, H Pr H 

We retain the assumption that the Prandtl number is unity, 
note that Re<— pV<r/n r) and substitute equation (25) into 
equation (42) to obtain 

F ’ < 43 > 

Now it is well known that at the high temperatures of 

interest here, the coefficient of viscosity p T varies nearly as 

the square root of the absolute temperature and is given by 
the relation 

Mr=2.31X10- 8 7V* 

If this expression is combined with equation (25) (neglect- 
ing T*), equation (43) may then be written 9 

^=6.8X10-^ V* (44) 


• The assumption of constant density certainly may invalidate this analysis for any quanti- 
tative study of the relatively “cold-wall” flows of Interest here. For the purpose of studying 
relative heat transfer it should, however, prove adequate, 

* The constant In equation (44) is obtained with the assumption of Incompressible flow In 
the stagnation region. The effects of compressibility and dissociation of the molecules of 
air in the region tend to increase the value of the constant by as much as a factor 2 In the 
speed range of interest in this report. For the comparative purpose of this report It is un- 
necessary to take these effects Into account. 


which, when combined with equations (2) and (13), becomes 


^ ^ (^PPy A 

^*=6.8X10-»y^ Vs 3 e 2 e Vmsin, s‘ (45) 


The maximum value of dHJdt can readily be shown to be 

■ V 

which occurs at the altitude 


0m sin e B Tr , ^ 

SfVJ 1 ' (46) 



corresponding to the velocity 

F,= F,r*a<0.85V* (48) 


The manner in which the heat-transfer rate per unit area at 
the stagnation point varies with incremental change in alti- 
tude from y z can be shown to be 


{dHJdt) it _ 
(dHJdt) 3 e 


2 e =F"'(0Ay) 


(49) 


The dependence of F nf (0A y) on $A y is shown in figure 4. 

Equation (46) applies only if y z is above sea level. If y 3l 
from equation (47), is negative, then the maximum heat- 
transfer rate per unit area at the stagnation point occurs at 
sea level and is 


dH^ / dH $ \ 

v dt / max \ dt )q 


6.8X10- 


6 jpo 


IV e 


Z(?[yp 0 A. 
S^rngin 9 E 


(50) 


DISCUSSION 

MOTION 

The motion study shows some important features about 
the high-speed descent of missiles through the atmosphere. 
The major assumptions of this analysis were that the drag 
coefficient was constant and the density varied exponentially 
with altitude. It was found that the deceleration due to 
drag was generally large compared to the acceleration of 
gravity and, consequently, that the acceleration of gravity 
could be neglected in the differential equations of motion. 
The flight path was then seen to be a straight line, the missile 
maintaining the flight-path angle it had at entry to the 
atmosphere. 

For most missiles, the maximum deceleration will occur 
at altitude. One of the most interesting features of the 
flight of such a missile is that the maximum deceleration is 
independent of physical characteristics (such as mass, size, 
and drag coefficient of the missile), being dependent only on 
the entry speed and flight-path angle (see eq. (17)). The 
missile speed at maximum deceleration (eq. (16)) bears a 
fixed relation to the entrance speed (61 percent of entrance 
speed), while the corresponding altitude (eq. (15)) depends 
on the physical characteristics and the flight-path angle but 
not on tl\e entrance speed. It is also notable that for a 
given incremental change in altitude from the altitude for 
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maximum deceleration, the deceleration and speed bear 
fixed ratios to the maximum deceleration and the entry 
speed, respectively (see fig. 4 and eqs. (19) and (18)), hence, 
the deceleration and speed variation with altitude can readily 
be determined. 

If the missile is very heavy, the calculated altitude for 
maximum deceleration (eq. (15)) may be fictitious (i. e., 
this altitude is negative) so the maximum deceleration in 
flight, which occurs just before impact at sea level, is less 
than that calculated by equation (17) and is dependent on 
the body characteristics as well as the entry speed and flight- 
path angle (see eq. (21)). However, the variation of speed 
and deceleration with altitude from the fictitious altitude 
given by equation (15) can still be obtained from figure 4. 

HEATING 

Total heat input. — In the heating analysis, a number of 
simplifying assumptions were made which should limit its 
applicability to the determination of relative values of heat- 
ing at hypersonic speeds. It is in this relative sense that the 
following discussion pertains. 

In considering the total heat transferred by convection 
to a missile, it is evident from equation (33) that the course 
the designer should take to obtain the least heating is 
affected by the value of the factor 

CppA =b (51) 

0m sin 0* 

To illustrate, first consider the case of a “relatively 
heavy” missile for which this factor is small compared to 
unity (the term “relatively heavy” is used to denote that 
the denominator involving the mass is very large as compared 
to the numerator involving the drag per unit dynamic 
pressure, C&A ) . Then 


\— e 


is small compared to 1. If this function is expanded in 
series and only the leading term retained, equation (33) 
becomes 


0 _ C/Sp.Vf 

40$\nd B 


(52) 


For the relatively heavy missile, then, the least heat will be 
transferred when C/S is a minimum — that is to say, when 
the total shear force acting on the body is a minimum. 
This result is as would be expected, if one notes that requir- 
ing #<<0 is tantamount to requiring the missile to be so 
heavy that it is retarded only slightly by aerodynamic drag 
in its motion through the atmosphere. Hence, the heat 
input to the missile is simply proportional to the shear force. 

Now Jet us consider the case when 1, or, in other 

words, when this missile is “relatively light.” In this event, 

C £fp 0 A. 

\— e < Smsin ^ ] 

and equation (31) can be approximated 


For the relatively light missile, then, the least convective 
heating is obtained when C/S/ C&A is a minimum. This is 
at first glance a rather surprising result, for it indicates 
that the heating is reduced by increasing the total drag, 
provided the equivalent frictional drag is not increased 
proportionately as fast. Physically, this anomaly is 
resolved if the problem is viewed in the following way: 
The missile entering the atmosphere has the kinetic energy 
imV K 2 but, if 


— CpP ° ' 4 - / y \ 2 

- 0m sin 9 S t \ 

* \vj 

is small, then nearly all its entrance kinetic energy is lost, 
due to the action of aerodynamic forces, and must appear 
as heating of both the atmosphere and the missile. The 
fraction of the total heat which is given to the missile is, 10 
from equation (33), 


1 fC/S\ 

2 \C d A) 

Thus, by keeping this ratio a minimum, as much as possible 
of the energy is given to the atmosphere and the missile 
heating is therefore least. 

In order to illustrate these considerations in greater detail, 
calculations have been made using the previously developed 
equations to determine the heat transferred by convection 
to a series of conical missiles. Two classes of missiles have 
been considered. Missiles in the first class were required 
to have a base area of 10 square feet. Missiles in the 
second class were required to have a volume of 16 cubic 
feet. Gross weights of 0, 1,000, 5,000, 10,000, and ® 
pounds have been assumed, and the entrance angle, 9 S , 
has been taken as 30° of arc in all cases. Missile heating, 
up to the tinre of impact, has been calculated as a function 
of cone angle for entrance speeds of 10,000, 20,000, and 30,000 
feet per second. In these calculations the pressure drag 
coefficient was taken as constant for a particular cone at the 
value corresponding to the entrance Mach number (a value 
of r*=500 o R was assumed throughout). These coeffi- 
cients were determined from reference 6 for cone angles of 
10° and greater. For cone angles less than 10°, reference 7 
was employed to determine these coefficients (base drag 
was neglected in all cases). The total drag coefficient was 
taken as the sum of the pressure drag coefficient plus the 
skin-friction coefficient, the latter coefficient being taken at 
its value for maximum total heat-input rate with altitude. 
The boundary layer was assumed to be wholly turbulent 
since the Reynolds number, based on length of run along 
the surface of a cone and local conditions just outside the 
boundary layer, was always greater than about 6X10 8 and, 
in fact, was of the order of billions for the more slender cones. 
Turbulent-boundary-layer data were obtained from refer- 
ences 8 and* 9, and Sutherland's law for the variation of 
viscosity with temperature was used in obtaining “equiv- 
alent flat-plate” heat-transfer coefficients. 


» Note that even that If all drag Is Motional drag, only half rhe heat is transferred to 
the body. The other half is contained in the boundary layer and is lert in the air in the 
body wako. 



Heal transferred, 0, foot-pounds 
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(a) V s = 10,000 ft/sec 

(b) F fi = 20,000 ft/sec 

(c) V E = 30,000 ft/sec 

Figure 5. — Convective heat transferred at impact to conical missiles 
of same area entering the earth’s atmosphere at an angle of 30° to 
the horizontal and velocities of 10,000, 20,000, and 30,000 ft/sec 
(base area=10 sq ft). 


Missile weight = <e lbs 





(c) Cone angle, degrees 


(a) Vg— 10,000 ft/sec 

(b) Vg= 20, 000 ft/sec 

(c) F* — 30,000 ft/sec 

Figure 6. — Convective heat transferred at impact to conical missiles of 
same volume entering the earth's atmosphere at an angle of 30° to 
the horizontal and velocities of 10,000, 20,000, and 30,000 ft/sec 
(volume= 16.34 cu ft). 

Missile heating calculated in this manner for the fixed-base- 
area and fixed -volume cones is presented in figures 5 and 6, 
respectively. Curves for missiles having densities greater 
than steel are considered improbable and are shown as 
dashed lines. It is clear that for both classes of bodies, when 
the missile is relatively heavy, the optimum solution is 
obtained by making C/S as small as possible (small cone 
angle easel and this optimum is accentuated with increase in 
speed. On the other hand, when the missile is relatively 
light, reduced heating is obtained by making C/S/C D A as 
small as possible (the large cone angle case). It is noted 
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(a) V g — 10,000 ft/sec 

(b) V g = 20,000 ft/sec 

(c) V M = 30,000 ft/sec 


Figure 8. — Maximum average rate of convective heat transfer to 
conical missiles of the same volume entering the earth’s atmosphere 
at an angle of 30° to the horizontal and velocities of 10,000, 20,000, 
and 30,000 ft/sec (voIume= 16.34 cu ft). 

temperature gradients through the shell that are intolerable 
even with the most highly conductive materials (copper, 
silver, etc.). 13 Thus some additional means of cooling, such 
as transpiration cooling, may, in any case, be required in 
this region. 

It was stated previously that pointed-nose bodies are un- 
desirable due, in part, to the fact that the local heat-transfer 
rate per unit area at the tip is excessive. The validity of this 
statement is demonstrated by the results of the analysis. 
It is clear (see eq. (44)) that since the local transfer rate varies 
inversely with the square root of the tip radius, not only 
should pointed bodies be avoided, but the rounded nose 


should have as large a radius as possible. The question 
then arises; if the nose radius is arbitrarily fixed, what course 
is available to the missile designer to minimize the problem of 
local heating at the stagnation point? From both equations 
(46) and (50), it is seen that for an arbitrary nose radius, if 
the mass, entry speed, and flight-path angle are fixed, then 
the only way to reduce the stagnation rate of heat input 
per unit area is to increase the product C&A. In fact, a 
relative stagnation-point heat-transfer rate per unit area, ^ 
may be expressed in terms of B (see eq. (51)), if it is defined 
as the ratio of the maximum stagnation-point heat-transfer 
rate per unit area for a given missile to the maximum rate 
the same missile would experience if it were infinitely heavy. 
For the infinitely heavy missile, the maximum rate occurs 
at sea level and is (see eq. (50)) 

6.8X10 TV 


so that from equation (50) 

3C i)P„-A 

^=e _29m5in9 *=e- »° (56) 

if the given missile also attains its maximum rate at sea level 
(i. e., y 3 =0; eq. (47)) ; whereas 


V ffm sin d_B __ 1 

3eC DPo A ~V3 '7B 


(57) 


if the given missile attains its maximum rate above sea level 
(eq. (46), y 3 positive). The variation of i with 1/5 is shown 
in figure 9. Clearly, the high pressure drag shape has the 
advantage over the slender shape in this respect. 

In order to illustrate these considerations in greater detail, 
again consider the constant-base-area and constant-volume 
cones discussed earlier. Assume the pointed tips of all the 



» See reference 1 for further discussion. 



12 


REPORT 1381— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


cones are replaced by spherical tips of the same radius c\ 
The relative effect of varying the cone angle on the stagna- 
tion-point heating can than be assessed by determining the 
variation of the product 



This product has been calculated for the various cones, as- 
suming C D to be unaffected by the addition of the hemi- 
spherical tip (the tip radius may be arbitrarily small), and 
the results are shown in figures 10 and 11. It is seen again 
that the missiles having large cone angle (high drag co- 
efficient) are considerably superior. 



(a) V b— 10,000 ft/sec 

(b) = 20,000 ft/sec 

( c ) Vb = 30,000 ft/scc 

Figure 10. — Maximum rate of convective heat transfer to the stagna- 
tion point of spherically tipped cones of the same base a^ea entering 
the earth's atmosphere at an angle of 30° to the horizontal and 
velocities of 10,000, 20,000, and 30,000 ft/scc (base area=10 sq ft). 





(a) V e — 10,000 ft/sec 

(b) V* = 20,000 ft/sec 

( c ) Vb = 30,000 ft/sec 

Figure 11. — Maximum rate of convective heat transfer to the stagna- 
tion point of spherically tipped cones of the same volume entering 
the earth's atmosphere at an angle of 30° to the horizontal and 
velocities of 10,000, 20,000, and 30,000 ft/sec (vo!ume = 16.34 cu ft). 

DESIGN CONSIDERATIONS AND CONCLUDING REMARKS 

In the foregoing analysis and discussion, two aspects of the 
heating problem for missiles entering the atmosphere were 
treated. The first concerned the total heat absorbed by the 
missile and was related to the coolant required to prevent 
its disintegration. It was found that if a missile were 
relatively light, the least required weight of coolant (and 
hence of missile) is obtained with a shape having a high 
pressure drag coefficient, that is to say, a blunt shape. On 
the other hand, it was found that if the missile were relatively 
heavy the least required weight of coolant, and hence of 
missile, is obtained with a shape having a low skin-friction 
drag coefficient, that is to say, a long slender shape. 

The second aspect of the heating problem treated was 
concerned with the rate. of heat input, particularly with 
regard to thermal shell stresses resulting therefrom. It was 
seen that the maximum average heat-input rate and, hence, 
maximum average thermal stress could be decreased by 



MOTION AND HEATING OF BALLISTIC MISSILES 


13 


using either a blunt or a slender missile, while missiles of 
intermediate slenderness were definitely to be avoided in 
this connection. The region of highest local heat-transfer 
rate and, hence, probably greatest thermal stress was 
reasoned to be located at the forward tip of the missile in 
most cases. This was assumed to be the case and it was 
found that the magnitude of this stress was reduced by 
employing a shape having the largest permissible tip radius 
and over-all drag coefficient; that is to say, the blunt, high 
drag shape always appears to have the advantage in this 
respect. 

These results provide us with rather crude, but useful, 
bases for determining shapes of missiles entering the atmos- 
phere which have minimized heat-transfer problems. If 
the over-all design considerations of payload, booster, et al, 
dictate that the re-entry missile be relatively heavy in the 
sense of this report, then it may be most desirable to make 
this missile long and slender, especially if the entry speed is 
very high (say 20,000 ft/sec or greater). Perhaps the slender 
conical shape is appropriate for such a missile. It seems 
clear, too, that the tip of this missile should be given the 
largest practicable nose radius in order to minimize the 
maximum local heat-transfer rate and hence maximum local 
shell stress problem. Even then it may be necessary to 
employ additional means to minimize the heat-transfer rate 
and, hence, thermal stress encountered in this region (e. g., 
by transpiration cooling). 

Let us now consider the case where the over-all design 
conditions dictate that the re-entry missile be relatively 
light in the sense of this report. This case will be the more 
usual one and, therefore, will be treated at greater length. 

A shape which should warrant attention for such missile 
application is the sphere, for it has the following advantages: 

1. It is a high drag shape and the frictional drag is only 

a few percent of the total drag. 

2. It has the maximum volume for a given surface area. 

3. The continuously curved surface is inherently stiff 

and strong. 

4. The large stagnation-point radius significantly assists 

in reducing the maximum thermal stress in the 
shell. 

5. Aerodynamic forces are not sensitive to attitude and, 

hence, a sphere may need no stabilizing surfaces. 

6. Because of this insensitivity to attitude, a sphere 

may purposely be rotated slowly, and perhaps even 
randomly 14 during flight, in order to subject all 
surface elements to about the same amount of 
heating and thereby approach uniform shell 
heating. 


u Note that if rotation is permitted, slow, random motion may be required in order to 
prevent Magnus forces from causing deviation of the flight path from the target. It should 
also be noted that at subsonic and low supersonic speeds gun-fired spheres, presumably not 
rotating, have shown rather large lateral motions in flight (see ref. 10). It is not known 
whether such behavior occurs at high supersonic speeds. 


On the other hand, the sphere, in common with other very 
high drag shapes may be unacceptable if: 

1. The low terminal speed cannot be permitted (e. g., 

because of excessive wind drift). 

2. The magnitude of the maximum deceleration is 

greater than can be allowed. 

The first of these disadvantages of the sphere might be 
minimized by protruding a flow-separation-inducing spike 
from the front of the sphere to reduce the drag coefficient 
to roughly half (see ref. 11). Stabilization would now be 
required but only to the extent required to counterbalance 
the moment produced by the spike. Special provision would 
have to be made for cooling the spike. 

Both of the disadvantages of very high drag shapes may 
however be alleviated by using variable geometry arrange- 
ments. For example, an arrangement which suggests itself 
is a round-nosed shape with conical afterbody of low apex 
angle employing an extensible skirt at the base. With the 
skirt flared, the advantages of high drag are obtained during 
the entry phase of flight where the aerodynamic heating is 
intense. Later, the skirt flare may be decreased to vary the 
drag so as to produce the desired deceleration and speed 
history. If the deceleration is specified in the equation of 
motion (see motion analysis), the required variation of drag 
toefficient with altitude can be calculated. 

The examples considered, of course, are included only to 
demonstrate some of the means the designer has at hand to 
control and diminish the aerodynamic heating problem. 
For simplicity, this problem has been treated, for the most 
part, in a relative rather than absolute fashion. In any 
final design, there is, clearly, no substitute for step-by-step 
or other more accurate calculation of both the motion and 
aerodynamic heating of a missile. 

Even from a qualitative point of view, a further word of 
caution must be given concerning the analysis of this paper. 
In particular, throughout, we have neglected effects of 
gaseous imperfections (such as dissociation) and shock-wave 
boundary-layer interaction on convective heat transfer to 
a missile, and of radiative heat transfer to or from the 
missile. One would not anticipate that these phenomena 
would significantly alter the conclusions readied on the 
relative merits of slender and blunt shapes from the stand- 
point of heat transfer at entrance speeds at least up to about 
10,000 feet per second. It cannot tacitly be assumed, 
however, that this will be the case at higher entrance speeds 
(see Appendix B). Accurate conclusions regarding the 
dependence of heat transfer on shape for missiles entering 
the atmosphere at extremely high supersonic speeds must 
await the availability of more reliable data on tiu* static and 
dynamic properties of air at the high temperatures and 
pressures that will be encountered. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. 28, 1953 
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APPENDIX A 


SYMBOLS 


D 

f 


referenco area for drag evaluation, ft 2 
body factor, dimensionless 
(See eq. (51).) 

drag coefficient, dimensionless 
skin-friction coefficient based on conditions just 
outside the boundary layer, dimensionless 
equivalent skin-friction coefficient, dimensionless 

(See eq. (28).) 

ft-lb 

specific heat at constant pressure, s j ug ojj 

i ft-lb 

specific heat at constant volume, s j ug 

F f ,F" ,F"' functions of 0A y, dimensionless 
(See eqs. (18), (19), and (49).) 
acceleration due to force of gravity 

( taken as 32.2 — 2 ) 
\ sec / 

„ . 4 ft-lb 

convective heat-transfer coefficient, jp gec 


c; 


c. 


Q 


h 

H 

k 

m 

M 

Nu 

Pr 

Q 

Be 

S 

T 

t 

V 


, . ft-lb 

heat transferred per unit area, 

, .. ft-lb 

thermal conductivity, sec (°R/ft) 

mass, slugs 

Mach number, dimensionless 
Nusselt number, dimensionless 
Prandtl number, dimensionless 
total heat transferred, ft-lb 
Reynolds number, dimensionless 
surface area, ft 2 

temperature (ambient temperature of air at 
altitude y unless otherwise specified), °R 
time, sec 

velocit y- 


x,y 

z 

& 


0 

1 

2 

3 

E 

l 

r 

s 

w 


horizontal and vertical distance from impact 
point, ft 

ft 2 

variable of integration, 

constant in density — altitude relation, ft 
(See eq. (2).) 

ratio of specific heat at constant pressure to 
specific heat at constant volume, 0 P lC v 
dimensionless 

increment 

distance within the shell measured normal to 
shell surface, ft 

angle of flight path with respect to horizontal, 
deg 

_ . . slugs 

coefficient of absolute viscosity, f tgec 

. J • , slug 

air density, 

radius, ft 

relative heat-transfer factor, dimensionless 
(See eqs. (56)' and (57).) 

SUBSCRIPTS 


conditions at sea level (y— 0) 
conditions at altitude t/i (eq. (15)) 
conditions at altitude y<i (eq- (38)) 
conditions at altitude y$ (eq. (4/)) 
conditions at entrance to earth’s atmosphere 
local conditions 
recovery conditions 
stagnation conditions 
wall conditions 

conditions within the shell of the missile 
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APPENDIX B 

SIMPLIFYING ASSUMPTIONS IN THE CALCULATION OF AERODYNAMIC HEATING 


As noted in the main body of the report, the heating 
analysis is simplified by making the following assumptions: 

1. Convective heat transfer is of foremost importance; 

.that is, radiative effects may be neglected. 

2. Effects of gaseous imperfections, in particular dis- 

sociation, may be neglected. 

3. Effects of shock-wave boundary-layer interaction 

may be neglected. 

4. Reynolds’ analogy is applicable. 

5. Prandtl number is unity. 

The restrictions imposed by these assumptions will now be 
considered in some detail. 

In assumption 1, two simplifications are involved; namely, 

(1) radiation from the surface of the body is neglected, and 

(2) radiation to the body from the high-temperature dis- 
turbed air between the shock wave and the surface is neg- 
lected. The first simplification may be justified on the 
premise that the maximum allowable surface temperature 
will be about the same for one body as compared with 
another, irrespective of shape, and, consequently, radiation 
away from the surface will be approximately the same. 
Hence, neglecting this form of heat transfer should not 
appreciably change the relative heating which is of principal 
interest in this paper. 

The second simplification of ignoring radiative heat trans- 
fer from the disturbed air to the body is not so easily treated. 
At ordinary flight speeds this form, of heat transfer is neg- 
ligible since it is well established that at temperatures not 
too different from ambient temperature, air is both a poor 
radiator and a poor absorber. At the flight speeds of 
interest, temperatures in the tens of thousands of degrees 
Fahrenheit may be easily obtained in the disburbed air flow, 
especially about the heavier blunt bodies. At these temper- 
atures it does not follow, a priori, that air is a poor radiator. 
Data on the properties of air at these temperatures are 
indeed meager. Hence, it is clear that calculations of 
radiative heat transfer from air under these conditions must, 
at best, be qualitative. Nevertheless, several such calcula- 
tions have been made, assuming for lack of better informa- 
tion that air behaves as a grey body radiator and that 
Wein’s law may be used to relate the wave length at which 
the maximum amount of radiation is emitted to the temper- 
ature of the air (this assumption, in effect, enables low- 
temperature data on the emissivity of air to be used in 
calculating radiation at high temperatures). In these 
calculations effects of dissociation in reducing the temper- 
ature of the disturbed air have also been neglected and 
hence from this standpoint, at least, conservative (i. e., too 
high) estimates of radiative heat transfer should evolve. 
The results of these calculations indicate the following: (1) 


Radiative heat transfer from the disturbed air to the body 
is of negligible importance compared to convective heat 
transfer at entrance speeds in the neighborhood of, or less 
than, 10,000 feet per second; (2) Radiative heat transfer, in 
the case of relatively massive blunt bodies, may have to 
be considered in heat-transfer calculations at entrance speeds 
in the neighborhood of 20,000 feet per second ; (3) Radiative 
heat transfer, in the case of relatively massive blunt bodies, 
may be of comparable importance to convective heat transfer 
at entrance speeds in the neighborhood of 30,000 feet per 
second. From these results, we conclude, then, that the 
neglect of radiative heat transfer from the disturbed air 
to the body is probably permissible for all except, perhaps, 
very blunt and heavy shapes at entrance speeds up to 20,000 
feet per second. However, this simplification may not be 
permissible, especially in the case of heavy blunt bodies 
entering the atmosphere at speeds in the neighborhood of, or 
greater than, 30,000 feet per second. 

In assumption 2, the neglect of effects of gaseous imper- 
fections, particularly dissociation, on convective heat transfer 
would appear to be permissible at entrance speeds up to 
and in the neighborhood of 10,000 feet per second, since at 
such speeds the temperatures of the disturbed air are not 
high enough for these imperfections to become significantly 
manifest. On the other hand, as the entrance speeds ap- 
proach 20,000 feet per second, temperatures of the disturbed 
air may easily exceed 10,000° Rankine, in which case appre- 
ciable dissociation may be anticipated, inside the boundary 
layer for all bodies, and inside and outside the boundary 
layer in the case of blunt bodies. The magnitude of these 
effects is at present in some doubt (see, e. g., the results 
of refs. 12 and 13.) Hence, for the present, the neglect of 
effects of gaseous imperfections on convective heat transfer 
is not demonstrably permissible at entrance speeds in the 
neighborhood of 20,000 feet per second or greater. 

In assumption 3, it has been shown by Lees and Probstein 
(ref. 14), and more recently by Li and Nagamatsu (ref. 15), 
that shock-wave boundary-layer interaction may signifi- 
cantly increase laminar skin-friction coefficients on a flat 
plate at zero incidence and Mach numbers in excess of about 
10. Lees and Probstein found somewhat the opposite effect 
on heat-transfer rate in the case of weak interaction. It is 
not now known how this phenomenon depends upon body 
shape or type of boundary layer. However, it is reasonable 
to anticipate that there will be some effect, and certainly if 
the skin-friction coefficient is increased in order of magnitude 
at Mach numbers approaching 20, as indicated by the 
results of Li and Nagamatsu for strong interaction, then 
the phenomenon cannot be presumed negligible. Hence, 
we conclude that from this standpoint, also, the convective 
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heat-transfer calculations of this report may be in error at 
entrance speeds of the order of 20,000 feet per second or 
greater. 

The assumption that Reynolds* analogy may be used to 
relate skin-friction and heat-transfer coefficient does not, 
especially in the light of recent work by Rubesin (ref. 16), 
seem out of line with the purposes of this paper, at least at 
entrance speeds up to and in the neighborhood of 10,000 
feet per second. However, it does not follow, a priori, that 
this assumption remains valid at substantially higher 
entrance speeds, especially in view of the imperfect gas and 
shock-wave boundary-layer-interaction effects already dis- 
cussed. 

The assumption of Prandtl number equal to unity would 
also appear permissible for the analysis of relative heating 
of missiles at the lower entrance speeds considered here. 
However, in view of the questionable effect (see again refs. 
12 and 13) of dissociation on Prandtl number, it is not clear 
that this assumption is strictly valid at the intermediate 
and higher entrance speeds treated in this report. 

From these considerations it is concluded that the simpli- 
fying assumptions made in the main heat-transfer analysis 
of this paper will not significantly influence the results at 
entrance speeds in the neighborhood of or less than 10,000 
feet per second. However, at entrance speeds in the 
neighborhood of and greater than 20,000 feet per second, 
these results must be viewed with skepticism. More accu- 
rate calculations of heat transfer at these speeds must, 
among other things, await more accurate determinations of 
both the static and dynamic properties of air under these 
circumstances. 
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